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Abstract 

We calculate the production of the gravitational waves from a double inflation 
model with lattice simulations. Between the two inflationary stages, gravitational 
waves with a characteristic frequency are produced by fluctuations of the scalar 
fields enhanced through parametric resonance. The wavelength of the produced 
gravitational waves gets extra redshift during the second inflationary stage and it 
can be in the observable range for the direct gravitational wave detectors. It is 
found that there is a possibility for the produced gravitational waves to be detected 
in the planned experiments. 
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1 Introduction 

The general relativity proposed by Einstein predicts the existence of the gravitational 
waves, which are the distortions of space-time geometry propagating through space as 
waves. Though it is indirect, there is strong evidence for the existence of the gravitational 
waves which is provided by the observation of the binary pulsar 1913 + 16 [1]. Yet 
up to date no direct detection has been made. In the world, there are lots of ongoing 
and planned experiments to detect the gravitational waves directly. Laser Interferometer 
Gravitational- Wave Observer (LIGO) [2] is working on the ground in USA. In Europe, 
Virgo [3] commenced operations in 2007 and Einstein Telescope (ET) [4] is planned. In 
Japan, KAGRA [5] is now under construction. It is also planned to construct the space- 
borne interferometers such as Laser Interferometer Space Antenna (LISA) [6], Deci- 
hertz Interferometer Gravitational Wave Observatory (DECIGO) [7, 8, 9], and Big Bang 
Observer (BBO) [10]. 

On the other hand, it is presumed that in the very early universe the vacuum energy 
dominates the universe and causes quasi-exponential expansion of the universe called 
inflation. Inflation explains the homogeneity and flatness of the universe. Wilkinson 
Microwave Anisotropy Probe (WMAP) [11] has measured anisotropics of the cosmic mi- 
crowave background (CMB) and successfully determined the cosmological parameters 
with high precision, which supports the prediction of the inflationary theories. However, 
it is impossible to prove directly the era before photons decouple from electrons, if we use 
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the result of the observation of cosmic microwave background only. In order to study the 
early universe before photon decoupling, we use the gravitational wave since the grav- 
itational wave interacts with other particles very weakly and preserves the information 
about the early universe. 

In this paper, we discuss the gravitational wave production from inflation. There 
are several primordial origins of the gravitational waves: quantum fluctuations of the 
metric [12], reheating [13], domain walls [14, 15], cosmic strings [16], etc. In this paper, 
we focus on the reheating as the origin of gravitational waves. In general, the wave 
length of the gravitational waves produced from reheating is so short [13] that we cannot 
detect the gravitational waves even in the future gravitational wave experiments with 
high sensitivity. In this work, however, we show that the double inflation model, which 
was originally proposed in [17] to solve the initial value problem for the new inflation, 
produces gravitational waves with the wavelength long enough to be proved with the 
planned detectors. 

In the double inflation, there are two inflationary stages. At the first stage hybrid 
inflation occurs, followed by the second stage at which new inflation occurs. In the 
intermediate stage, the inflaton and waterfall fields oscillate around the minimum of the 
potential in the same way as the reheating. In this paper, we adopt the smooth hybrid 
inflation as the first inflation [18], in which the gauge symmetry of the waterfall fields is 
broken even during the inflation, thus no harmful topological defects are produced. In the 
previous study [19], it was found that in this model fluctuations of the scalar fields are 
enhanced by the parametric resonance and then primordial black holes are formed at a 
characteristic scale. In this paper, we pay attention to the production of the gravitational 
waves from these fluctuations. Since the parametric resonance is a nonlinear phenomena, 
we calculate the production of the gravitational waves with lattice simulations. 

The organization of this paper is as follows. We review the double inflation model in 
Sec. 2. In Sec. 3, we calculate the amount of the produced gravitational waves by using 
analytical consideration and lattice simulations. We estimate the present density of the 
gravitational waves and the peak frequency in Sec. 4. Finally, Sec. 5 is devoted to the 
conclusion. 

2 Smooth hybrid new inflation 

In this section, we review the smooth hybrid new inflation model in supergravity proposed 
in [18]. This model has two inflationary stages, where the smooth hybrid inflation [20] 
takes place first and then the new inflation [21] follows. We assume that the e-fold number 
of the new inflation is smaller than 60. Thus, the density fluctuations on large scales are 
produced during the smooth hybrid inflation. In this section, we set the reduced planck 
mass Mpi ~ 2.4 x 10 18 GeV equal to unity unless otherwise stated. 
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2.1 Smooth hybrid inflation 

First, we discuss the smooth hybrid inflation model. This model has three superfields. 
One is the inflaton field S and the others are waterfall fields \l/ and ^f. This model is 
based on U(1)r symmetry. The superpotential is given by 

w 

S has R charge 2 and fifty) has R charge 0. The i?-invariant Kahler potential is given by 

K H = \S\ 2 + \ty\ 2 + \ty\ 2 . (2) 

We also introduce the U(l) gauge symmetry under which the waterfall fields are trans- 
formed as ty — > e tS ty and ty — > e~ l5 ty. Using the U(l) phase rotation and D-term flat 
condition, we take ty = ty hereafter. 

From the superpotential and the Kahler potential, we can write down the scalar po- 
tential. Using the R symmetry, we can change the complex scalar field S into real scalar 
field, a = V2ReS. We assume that only one degree of freedom of ty dominantly describes 
the inflation dynamics. Hence we change the complex scalar field \1/ into real scalar field, 
ip = 2Rety. Neglecting higher-order terms, we can write down the scalar potential for 
a < 1 

where M is a cut-off scale, a is the inflaton and ip is the waterfall field. During inflation, 
ip is located at the local minimum of the potential, 

if for a » v^M, (4) 
2V/xM for a < y/flM. 

Since the gauge symmetry of the waterfall fields ty and ty is already broken during inflation, 
no topological defects are formed in this model. When a is larger than y r jlM, the effective 
potential is given by 

4 I", 2 fi 2 M 2 , 



V H (a) ~ fi 



27 a 4 + 8 



(5) 



The third term reflects the effect of the supergravity. Since the effective potential is 
dominated by the false vacuum energy /x 4 , the Hubble parameter is given by 



To consider the dynamics of the inflation, we differentiate the effective potential as 

(7) 



dV H 4 
^7 = /i 



8 (/iM) 2 + 1^ 



27 a 5 2 
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The second term, which reflects the effect of the supergravity, dominates the inflation 
dynamics when inflaton a is larger than ad given by 

<r d =(^y\»Mf\ ( 8) 

Let us define cr, as the value of the inflaton field when pivot scale ko = 0.002 [Mpc -1 ] 
leaves the horizon. The slow-roll condition is broken when \r)\ = ||^| ~ 1 and inflation 
ends at a ~ a r where 
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1/6 



ac= {^) (/iM)1/3 ' (9) 
Then, the number of e-folds Nh from cij to a c is written as 

N R = / da— = / da — — -5 . (10) 

27 (j ' 2 

The amplitude of the curvature perturbation in comoving gauge 71 and the spectral index 
n s at the pivot scale are written as 

WH _L^^w = — 1 



8 {jMY A- Irr 3 
27 ct 5 2 « 



2nVS \V'(a t )\ 2n ^ 
and 

80 (uM) 2 „ . , 

n s -l = -- K -^L + 3al (12) 
27 a? 

In this model, when the value of the inflaton at the pivot scale <Ji is larger than er^, the 
spectral index becomes blue. 

2.2 New inflation 

Second, we describe the new inflation model [21] we adopt. This model has one superfield 
$. $ has R charge | and U(1)r symmetry is dynamically broken down to a discrete Z 8 r 
at a scale v. The superpotential is given by 

W N = v 2 <$> - |$ 5 , (13) 

where g is a coupling constant. The R- invariant Kahler potential is given by 

K N = \<t>\ 2 + ^\<S>\\ (14) 

Here, Cn is a constant of order 1. 

The scalar potential yields a vacuum, 



2\ 1/4 



($) ~ - . (15) 
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At this vacuum, the potential has the negative vacuum energy as 

2 /,2\l/2 



V((<Z>)) --3e KN \W N [(^}] 



4 \ 4 / r~ 

v 



~-3 1 ; > - ; 

5/ V 



(16) 



We assume that the negative vacuum energy is cancelled out by a supersymmetry-breaking 
effect which gives a positive contribution Ag USY to the vacuum energy. Thus, we have a 
relation between v and the gravitino mass m 3 / 2 as 



A susy 4 2 / v 



2\ 1/4 



(17) 



To discuss the dynamics of the new inflation, we change the complex scalar field $ 
into the real scalar field: = \/2Re&. Neglecting the higher-order terms, we obtain the 
scalar potential of the new inflation as 



V H [(f)] ~v - —v (j) - -v (j) + — < 



We differentiate the scalar potential as 



,2j,3 , fi 1 /7 



(h> =-C N v 2 <P-2gv'<p A + »-<P<. 



(18) 



(19) 



When is smaller than a certain value 0^, the first term dominantly controls the inflation 
dynamics. On the other hand, when is larger than 0^, the second term dominantly 
controls the inflation dynamics. Here, 0^ is given by 



^9 



(20) 



New inflation occurs when inflaton is smaller than C . C is the value at which slow-roll 



condition is broken, \rj\ 



1 1 V" 1 

- 

2 I Viv 



1, and is given by 



6g 



:i-c 



(21) 



Then, the number of e-folds iV ncw is written as 

N 

1 "new 

Here we assumed that Cn < \- 



10, 1-4C N 
in - — h 



C N <pi 2C N (1 - C N ) 



(22) 
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2.3 Oscillatory phase 

Here, let us consider the oscillatory phase from the end of the smooth hybrid inflation to 
the start of the new inflation. 

After the end of the smooth hybrid inflation, a and ip roll down toward cr m i n = and 
VVin — 2-y/yuM where the potential has the minimum. Around the potential minimum, 
their effective masses m a and are given by 



(23) 

^ + 16/i 4 > H. 

Thus, a and ip oscillate around their respective minima and the total energy density 
decreases as a~ 3 . Eventually, the false vacuum energy of the new inflation v 4 dominates the 
total energy density of the universe, and then new inflation starts. During this oscillatory 
phase, the effective masses of the inflaton and waterfall fields change periodically through 
their mutual couplings, which leads to the rapid amplification of the field fluctuations 
5a and dip called the parametric resonance [22, 23, 24]. In the present model, the most 
amplified fluctuations have wavenumber /c pea k [19] given by 

fc peak ~ 0.35m ff . (24) 

We expect that these enhanced fluctuations produce gravitational waves at the scale /c pea k- 
Next, we consider the initial condition for the new inflation. Here, interactions between 
the two inflation sectors become important. During the smooth hybrid inflation, the 
leading term of the interaction between a and <p is given by 

^ n t = ^V + /iVa0. (25) 

Thus, at the end of the smooth hybrid inflation, the value of the inflaton of the new 
inflation is given by 

v 2 

(pram = " c . (26) 

During the oscillatory phase, the potential energy of the scalar fields averaged over one 
oscillation time is the half of the total energy density of the universe, which leads to the 
effective mass of <p, 

ml = // = V = \E\ (27) 

Taking into account H = \H 2 1 and using eq. (27), one can find that the amplitude of (p 
decreases as a _3//4 . Thus, the initial condition for the new inflation is estimated as 

3 

0ini ^ -\(Tc- (28) 
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Taking into account this oscillatory phase, and assuming the ordinary thermal history 
after the reheating, we estimate the total number of e-folds iV tot a s 

iVtat = N H + N new = 57.6 + ln/i + i In (^) , (29) 

where we set the reheating temperature as Tr. 



2.4 Parameter search for Oi 

In this subsection, we search for the allowed region of <jj at each n s . At the pivot scale 
ko, the spectral index n s and the curvature perturbation 1Z are written in terms of <Tj, the 
energy scale of the smooth hybrid inflation fi and the cut-off scale M. Then, using the 
WMAP constraint 1Z ~ 4.7 x 1CT 5 , we can write down /i and M in terms of n s and <7j as 



H = 7.2 x 1(T J J 8 erf - (n a - 1) tr f , (30) 



A f =8 i^/ o 3 :?-'"--;> . (3D 

»' ° of - (n a -l)<7i 

The value of the inflaton a c , at which the smooth hybrid inflation ends, is also written in 
terms of cr, and n s as 

<7 C = 0.89 < T i [3<7 f 2 -(n a -l)] 1/fl . (32) 

The e-folding number of the smooth hybrid inflation Nn and that of the new inflation 
N new are written as 

A'u ■ / da- ^ 



8 (MM) 2 
27 ct 5 



(33) 



c?cr — — - — : , 



0.89 aj 3 ^-(n s -l)l 1/6 1 *f foHn.-p) CT 3 

10 - 5 



A^ncw = iVtot - N U 



54.3 + - In [8^ 3 - (n s - 1) a t ] 



dc —, — ? 

(1.89 CTi [3 CT ?-(n s -l)l 1/6 1 °f(3°f-(tt.-D) CT 3 
L J 10 <r 5 ^2 



(34) 



where we have used (29) and (30) for N tot and Thus, the parameters of the smooth 
hybrid new inflation are written in terms of <7; and n s . 

Since the effective potential in the supergravity is uncontrollable over the planck mass, 
(7j and fi must be smaller than the planck mass. In addition to this, <7j must be larger 
than a c . This condition guarantees that the hybrid inflation occurs. After the hybrid 
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inflation, the fluctuations of the scalar fields are amplified through the parametric res- 
onance, and then these fluctuations lead to a sharp peak at the scale corresponding to 
k~ eak in the spectrum of the curvature perturbation. The observations of the large scale 
structures have not seen such enhancement of the curvature perturbations so that the 
scale corresponding to k~^ ak should be sufficiently small. Thus, Nn must be larger than 
about 10. Futhermore, since we assume that the new inflation occurs after the hybrid 
inflation, iV new should be larger than 0. In addition to these conditions, the cut-off scale 
M must be larger than fi. These conditions are summarized as 

(i) o c < oi < 1 

(ii) n < 1 

< (35) 

(1) fi< M 

k (iv) N u > 10, N new > 

Substituting (30) — (34) into this condition, we can obtain the maximum and the minimum 
values of crj. 

Fig. 1 shows the maximum and the minimum values of Oi at each value of n s . Using 
(30) and (31), we can get the corresponding maximum and the minimum values of \x and 
M which are shown in Fig. 2. 



3 Production of the gravitational waves 

During the oscillatory phase, the fluctuations of the scalar fields grow exponentially 
through the parametric resonance. These growing fluctuations can be a source for grav- 
itational waves. In this section, we calculate the amount of the produced gravitational 
waves. First, we roughly calculate the energy density of the gravitational waves. Second, 
we introduce the formalism with which we calculate the energy density and the spectrum 
of the gravitational waves produced from the scalar fields. Finally, we perform the lattice 
simulations and calculate the evolution of the fluctuations of the scalar fields. Using the 
result of the lattice simulations, we estimate the energy density and the spectrum of the 
gravitational waves. 

3.1 Order estimate 

In this subsection, we roughly calculate the ratio fi gWiP of the energy density of the grav- 
itational waves to the total energy density based on the naive estimation given by Felder 
and Kofman [13, 25]. Here and hereafter the subscript p means the epoch, at which 
gravitational waves are produced. 

Through the parametric resonance, the fluctuations of the scalar fields with the typical 
length scale -R pea k — &p ea k are enhanced most. These fluctuations produce the gravita- 
tional waves. The enhanced fluctuations have self gravitational potential energy E gw . We 
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Figure 1: Maximum and Minimum values of <7j at each value of n s obtained from the 
condition (35). The upper panel shows the maximum and minimum values of <7j, satisfying 
all conditions of (35). The lower panel shows each maximum and minimum value of crj, 
satisfying each condition given by (35). The red line is the upper limit satisfying a c < <Ji. 
The green line is the upper limit satisfying N-& > 10. The purple line is the lower bound 
on CTj coming from [i < M. The lower limit on <jj from iV new > is represented by the sky 
blue line. The blue line shows the lower limit for existence of fi and M. 
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Figure 2: Maximum and Minimum values of \i and M at each value of n s obtained 
from the condition (35). The upper panel shows the maximum and minimum values of 
fi, satisfying all conditions of (35). The lower panel shows the maximum and minimum 
value of M, satisfying all conditions of (35). 
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can write down E gw at t p as 



£ gw ~ G {ptot f^ k)2 * p tot , P H*Rl eak , (36) 

-'kpeak 



where ptot,p is the total density and H p is the Hubble parameter. We approximate that 
the order of the energy of the gravitational waves is comparable to the self gravitational 
potential energy of the enhanced scalar fields E^. Under this approximation, the energy 
density of the gravitational waves p gw ,p is estimated as 

Pgw,p ~ -^3^ ~ Ptot,p#p-Rp eak - ( 37 ) 

Thus, the ratio of the energy density of the gravitational waves to the total energy density 
fi gWjP is given by 

^gw,p = ^ SW ' P ~ {RpeakHp) 2 . (38) 
Ptot,p 

In this paper we assume that the density of the produced gravitational waves is given by 
the above equation with introducing a numerical constant as 

fi g w,p = OL {RpeakHp) 2 , (39) 

where a will be determined by the numerical simulations presented in Sec. 3.3. For the 
smooth hybrid new inflation model, from (24), -R pea k is written as 



1 _ 1 M 
peak " 035^" 035 V8^' 



Using this relation, we can rewrite (39) as 

*W = 0.34 « ) ■ (41) 

3.2 Formalism for calculating gravitational wave spectrum 

Here, we review the formalism developed by Dufaux et al. [27] with which we calculate the 
spectrum of the gravitational waves produced from the fluctuations of the scalar fields. 

We consider the Friedmann-Robertoson- Walker metric including the linear perturba- 
tion in the spatial metric Ttij. 

ds 2 = a 2 (r) [-dr 2 + (<J y + ^(r, £)) dx L dx j ] , (42) 

where r is the conformal time defined as dt = adr. The amplitude of the gravitational 
waves hij is defined as the transverse-traceless part of %ij 

hij(T,x) = a(r)nj^(r,x). (43) 



11 



The superscript TT means the transverse-traceless part. The evolution of hy is described 
by the linearized Einstein equation as 



h'lfax) - V 2 ^(r,x) = lQnGa(r)T^ T (r,x 



(44) 



Assuming that the typical physical scale is much smaller than the horizon scale, we 
neglected the terms including a". Here, T? t (t, x) is computed by applying the projection 
operator in the momentum space as 



7f T (r, k) = A ijM (k)T tJ (r, k) = A ijM {<9 fc 0<9,0} (r, k), 



Pij\k) — Sij kikj, 



(45) 



where {<9fc0<9/0} (r, k) is the Fourier transform of dk4>(r, x)di<p(T, x) . 

We assume that the source term is nonzero only for the time interval Tj < r < tj. In 
this case we can solve (44) by using the Green's function and obtain 



hij(r, k) = Aij(k)sin [k (r — Tj)] + Bij(k)cos [k (r — T/)] for r > 77, 



(46) 



where 



Aij[k) 



-16ttG n 



dr' cos [Jfe ( Tf - r')] a(T')T^ T (T', k) 



k 



16nG n 



k 



(47) 



dr' sin [k (r/ — r')] a(r')T i ■ (r', k). 



— * 

Here, hij(r, k) is the Fourier transform of hij(r,x). 

The energy density of the gravitational waves is given by 



Pgw{r) 



327rGa 4 ( 



T 



(^.(r,f)^ 7 (r,f)) 



32irGa 4 (r)V J (2tt 



d 3 k 



KAT,k)hUT,k), (48) 



where (• • -)y denotes average over spatial volume 1/. Taking the time average over a 
period of the oscillations, we obtain p gw ( r ) as 



Pgw(r) 



4ttG 



\r)Vj (2tt) 



E 



dr' cos (fcr') a (r') T^ T (r', fc) 



+ 



dr' sin (kr) a (r') T^ T (r', fc) 



Thus, the spectrum of the gravitational waves just after the production is given by 

1 d Pgw (T f ,k) _ AG 2 



ptot dink 3irVa 4 (T f )H 2 (T f ) 



Sk(Tf), 



(49) 



(50) 
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where 



s k (r f ) = k 3 I dn k J2 



dr' cos (kr) a (r') T^ t (t', k) 



+ 



where dQ k = d cos 



dr' sin (kr) a (r') T? t (t', k) 



(51) 



3.3 Lattice simulation 

The parametric resonance is a nonlinear phenomenon. Thus, in order to estimate the 
spectrum and energy density of produced gravitational waves precisely, we need to perform 
the lattice simulations to calculate the time evolution of the fluctuations of the scalar fields. 

We perform the lattice simulations with three dimensional comoving box taking the 
cosmic expansion into account, where the scale factor is detrmined by the Friedmann 
equation. In order to resolve the oscillations of the inflaton, we set the time step dt as 



dt 



1 1 1 

10 rn„ m r 



(52) 



We set the step number of the time as 1500, where the ratio of the scale factor at the 
end of the simulation to that at the beginning is 3.1 l . By the end of the simulation, 
the production of the gravitational waves stops. The linear analysis [19] shows that the 
enhancement of the scalar fields occurs at the spacial scale -R pea k ~ 1/ (0.35m CT ). Taking 
this into account, we set the box size of the simulation V = L 3 and the grid number iV gri d 
as L = 100/m CT and A^ gri d = 128 3 . Then, L > -R pca k and L/N^ d < -R pea k are satisfied even 
if we consider the growth of the scale factor. Thus, we can resolve the peak scale -R pea k 
during the simulation. We choose the physical parameters as follows 



fi = 2.04 x 10 



-3 



5.0 x 10" 4 , g = 2.0x TO" 5 , and C N = 0.04, 



(53) 



which were used in [19] to analyze the oscillatory phase in the linear perturbation regime. 
We set the initial condition for the fluctuations of the scalar fields as random variables. 

The inflation model we consider here has three scalar fields a, ip and 0, where a is 
the inflaton of the hybrid inflation, if) is the waterfall field of the hybrid inflation and 
<f) is the inflaton of the new inflation. We perform the simulation using these three real 
scalar fields. 2 We set the initial conditions for the scalar fields as the value just after the 

1 If we continue the simulation for a larger time step, the energy density of the universe is dominated 
by the constant energy density v 4 and the scale factor increases exponentially with time, which determines 
the upper limit of the simulation time. 

2 More precisely, this model has three complex scalar fields S, $ and 4>. The real part of S, *f> and 
$ corresponds to a, ip an d 4>- The imaginary part of each complex scalar fields does not affect on the 
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Figure 3: Time evolution of the occupation number n/ from t = 30[l/m CT ] to t = 
150[l/m CT ] for every 20[l/m cr ] time interval. The top, middle and bottom panels show the 
time evolutions of n/ for a, ip and <f), respectively. 
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end of the first inflation including random fluctuations at each grid. Using the lattice 
simulations, we calculate the occupation number of the scalar fields in momentum space 
rif as 

n f = u k (m^ + |M 2 ) , ul = - 2 +m): f e {*, V, <P}, (54) 

V u k J a 

where fk is the Fourier component of /. Fig. 3 shows the time evolution of n/ for the 
scalar fields. We can see the enhancement of the scalar fields a and ip. This enhancement 
is induced by the parametric resonance. 

Using the result of the lattice simulations, we calculate the spectrum of the gravita- 
tional waves. Fig. 4 shows the time evolution of the spectrum of the gravitational waves 
produced from three real scalar fields. We can see that the spectrum of the gravitational 
waves produced just after the first inflation has a sharp peak at the scale of the inflaton 
mass 77V. This corresponds to the scale at which the fluctuations of the scalar fields 
grow through the parametric resonance. Furthermore, we perform the simulations for 30 
realizations with different initial random values of the scalar field fluctuations. Using the 
result of these 30 realizations and (50), we estimate the amount of the gravitational waves 

^gw I Lattice 

fi g w|Latticc = (2.0 ±0.3) x KT 3 . (55) 

With the lattice simulations, we can precisely calculate the energy density and the 
spectrum of the gravitational waves produced during the oscillatory phase. However, 
because of the cost of the simulation time, we cannot perform the lattice simulation for 
all parameters allowed in this model. In order to estimate the amount of the gravitational 
waves for various parameters, we use (41) determining the numerical constant a from the 
result of the lattice simulations. For the parameters given by (53), (41) leads to 

^gw (estimate ~ « X 7.8 X 10" 4 . (56) 

Comparing (55) with (56), we can determine the a as 

a ~ 2.6 ±0.4. (57) 
Thus, the amount of the gravitational waves just after the production is given by 

fi gw , p ^ (0.88 ±0.14)^. (58) 



dynamics of inflation directly. However, there is a possibility that the fluctuations of their imaginary 
parts are enhanced through the parametric resonance and the shape of the gravitational waves becomes 
different from that simulated only by three real scalar fields. We have performed the simulations with 
these three complex scalar fields, from which we confirmed that the spectrum is not altered by the 
fluctuations of the imaginary part. 
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GW from 3 real scalar field 




0.1 1 10 100 

k[2jt/L] 

Figure 4: The time evolution of the spectrum of the gravitational waves from t = 
30(1/7710-] to t — 150[l/m a ] for every 20[l/m CT ] time interval. Here, the vertical axis 
corresponds to Sk defined in (51). 



4 Present density of gravitational waves 

After the oscillatory phase, the universe is dominated by the false vacuum energy v A 
and the new inflation starts. After the new inflation, the universe experiences the re- 
heating followed by the standard thermal history. Because of the cosmic expansion, the 
gravitational waves are red-shifted and their energy density decreases. Considering these 
effects of the expansion, we calculate the present peak frequency of the gravitational wave 
spectrum / gWj o, and the amount of the gravitational waves f2 gw ,o- 



4.1 Peak frequency of gravitational waves 

As a result of the lattice simulation, it is found that the spectrum of the gravitational 
waves has a sharp peak at / gw ,p — & P eak/27r, 

/gw ' p= 27 035^ ' (59) 

This frequency decreases by the cosmic expansion and the present frequency / gW]0 is given 
by (cip/ao) f gw ,p, where a p is the scale factor at the production of the gravitational waves. 
We assume that the scale factor a p is approximately equal to the scale factor at the end 
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of the smooth hybrid inflation an- Then, the ratio a /a p is given by 



/ rp \ -1/3 / \ 4/3 

Thus, the present peak frequency of the gravitational waves / gWj o is given by 

A— *^ (61) 

4.2 Relic density of gravitational waves 

In this subsection, we consider the effect of the cosmic expansion on the energy density of 
the gravitational waves. Just after the production, the energy density of the gravitational 
waves decreases as a -4 , like radiation. After the wavelength becomes larger than the 
horizon during the new inflation (horizon exit), the amplitude of the gravitational waves 
stops decreasing, and the energy density decreases as a~ 2 . After the end of the new 
inflation, the gravitational waves re-enter the horizon (horizon entry), and the energy 
density decreases as a~ 4 again. Thus, the energy density of the gravitational waves today 
P g w,o is given by 

\ — 4 / \ 2 

\ / Gentry 



p g w,o = — Pgw, P - (62) 

V OH / V Oexit / 

Here, a ex it is the scale factor at the horizon exit. a cntry is the scale factor at the horizon 
entry. 

If iVnew is sufficiently small, the wavelength never exceeds the horizon during the new 
inflation. We call this case "No horizon exit" . As a { becomes larger, N Dew becomes larger 
and horizon exit occurs during the new inflation. After the new inflation, the horizon 
begins to grow, and the gravitational waves re-enter the horizon. If iV new is small, the 
horizon size at the end of the new inflation is not so large compared with the wavelength 
of the gravitational waves. Thus, the gravitational waves re-enter the horizon soon after 
the new inflation, i.e., before the reheating is completed. We call this case "Reheating 
time horizon entry". As N ncw becomes larger, the horizon entry occurs later, i.e. during 
the radiation dominated era. We call this case "Radiation time horizon entry". When 
iV ncw becomes much larger, the horizon entry occurs during matter dominated era. We 
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call this case "Matter time horizon entry". In each case, we can write down the fi gWi o, as 

-4/3 

( ) 



No exit /i 2 n gw , , no = 4.5 x 10~ 17 e" 4Af - (x^) 473 (^) ft 
Reh. entry /^ gw , , r ch= 6.4 x 10~ 19 e 2 ^™ (V) 4 (^) 4/3 f^)^ (ff) 3 ft gw , P 
Ra. entry /i 2 fi gw ,o, ra = 4.0 x 10~ 6 (j)* (f|)' ft gw , p 

Ma. entry ^,0,™= 8.7 x 10~ 55 e 2JV - (j)* (j^f* (^'^ (^) 3 *W 

(63) 

The epochs at which the horizon entry occurs are categorized into four cases under three 
conditions CI, C2 and C3. If CI is satisfied, the horizon entry occurs after the reheating. 
If C2 is satisfied, the horizon entry occurs after the beginning of the radiation dominated 
era. If C3 is satisfied, the horizon entry occurs after the matter-radiation equality time. 
These conditions are given by 



pi 



Reh. entry : CI e*— > 1.7 (j^j (j 
Ra. entry : C2 ^ > 1.6 x 10 6 (j^)' 1 (^)^ (fg)~ 1/2 (64) 
Ma. entry : C3 e»™ > 3.2 x 10 2 * (^) 1/3 (^f* 



4.3 fgwfi and n gW; o 

In this subsection, we calculate the relation between the amount of the gravitational waves 
today fi gw ,o and the peak frequency of the gravitational wave spectrum /gw,o- In (61) 
and (63), /gw,o an d ft gw ,o are written in terms of iV new , /1, v, M and Tr. Furthermore, 
in (34), (30) and (31), iV new , /1, M are written in terms of a^ and n s . Here, we assume 
that the energy scale of the new inflation v is approximately equal to that of the hybrid 
inflation /i, and the reheating temperature Tr is 10 9 GeV. For // ~ v the density of the 
gravitational waves is maximal so that our estimation gives the upper bound. With these 
assumptions, / gWj o and fi gw ,o are written only by Oi and n s as 

In /gw.o = -39.3 - ^ In a % - i In [3a 2 - (n s - 1)] + N H , (65) 

and 

' In fi gWi0 ,no = In a - 250 + In o { + § In [3<7< - (n s - 1)] - 2 In [8a 2 - (n, - 1)] + 4iV H 
lnfi gWi0ireh = In a + 70.1 + 12 In a* + 2 In [3a, - (n s - 1)] - 2iV H 
lnfi gWi o, ra = In a - 15.12 + 9 In a; + fin [3a 2 - (n s - 1)] 
In fl gw ma = In a - 12.5 + 12 In <n + 2 In [3a 2 - (n s - 1)] - 2jV H 

(66) 
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Figure 5: Relations between f2 gWj0 and / gWj0 for given spectral indeces n s = 
-0.04,-0.038,-0.03, and 0.1. 



where domain of is constrained as shown in Fig. 1. Thus, for a given spectral index n s 
we can obtain f2 gw> o as a function of the peak frequency / gw ,o- Fig- 5 shows this relation. 
At low peak frequencies / gw ,o ^ 40 Hz, the horizon entry occurs during the radiation 
dominated era. At high peak frequencies 40 Hz < / gWj0 < 10 5 Hz, the horizon entry 
occurs during the reheating era. At higher peak frequencies / gw ,o ^ 10 5 Hz, the horizon 
exit does not occur. In this model, as <jj becomes larger, \x and fi gWiP becomes larger. In 
addition to this, as <jj becomes larger, N new becomes larger. Hence, / gWi0 becomes smaller 
as shown in (61). When the wavelength is larger than the horizon, the amplitude of the 
gravitational waves does not decrease. Because of this effect, fi gWjP increases. Thus, the 
curves of the plot have negative slope. When horizon exit does not occur, the amplitude 
of the gravitational waves always decreases during the new inflation. This effect decreases 
the amount of the gravitational waves more efficiently than in the case where horizon exit 
occurs, and the slope of the plot becomes positive at higher frequencies. 

5 Conclusion 

In this paper, we have studied the production of the gravitational waves in the smooth 
hybrid new inflation model [18, 19]. After the smooth hybrid inflation, the inflaton and 
waterfall fields for the hybrid inflation start to oscillate and their fluctuations grow ex- 
ponentially through the parametric resonance, which leads to the efficient production of 
gravitational waves. We performed the lattice simulation and calculated the amount of 
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Figure 6: Sensitivity curves of planned detectors, ultimate DECIGO, BBO, LISA, ET, 
and KAGRA. The pink line shows the spectra of the gravitational waves in the present 
model with /i = 3.8 x 10" 4 , M = 4.1 x 1(T 2 and n s - 1 = -0.035. The blue line shows 
the spectrum of the gravitational waves from the single hybrid inflation. The dotted line 
shows the upper bound of the tensor to scalar ratio. 

the gravitational waves whose spectrum has a peak at the scale predicted by the linear 
analysis [19]. 

From the result of the lattice simulation, we found the relation between the present 
gravitational energy density Q gw> o and the peak frequency f gw . Fig. 6 shows the sen- 
sitivities of the planned detectors and gravitational wave spectrum predicted from in- 
flation models. This figure shows that the gravitational waves produced in this model 
can be detectable in the future experiments in some parameter space. For example, the 
pink line shows the spectrum of gravitational waves produced in a set of parameters 
[i = 3.8 x 10" 4 , M = 4.1 X 10" 2 , n s - 1 = -0.035, where n s is the best-fit value of 
the WMAP 7 yr result. In this parameter set, the ultimate DECIGO can detect the 
gravitational waves. 
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